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. In the marriage problem, a variant of the bi-parted matching problem, each member has a 'wish- 

■ list' expressing his/her preference for all possible partners; this list consists of random, positive real 

Q\ | numbers drawn from a certain distribution. One searches the lowest cost for the society, at the risk 

of breaking up pairs in the course of time. Minimization of a global cost function (Hamiltonian) is 
performed with statistical mechanics techniques at a finite fictitious temperature. 

The problem is generalized to include bachelors, needed in particular when the groups have 
different size, and polygamy. Exact solutions are found for the optimal solution (T = 0). The 
entropy is found to vanish quadratically in T. Also other evidence is found that the replica symmetric 
solution is exact, implying at most a polynomial degeneracy of the optimal solution. 

Whether bachelors occur or not, depends not only on their intrinsic qualities, or lack thereof, but 
£h ' also on global aspects of the chance for pair formation in society. 
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The stable marriage is a well-known optimization problem "women" , and pairs are formed between them. Many 
other applications exist, such as the job-market (employers and employees), the housing market, etc. Each man ranks 
the women according to his preference, number one being his most preferred candidate; likewise, each woman ranks 
the men. The cost function for a configuration of marriages is the sum of the rankings of the actual partners on the 
two lists. In the past most studies have been concerned with algorithms for finding stable pairing solutions. Indeed, 
given the rules according to which marriages may be ended and started at new, a stable situation may be reached 
where it brings no benefit to break up existing relations for starting other ones. For a recent study of this subject, 
_l | see Omero et al [|| . 

. Since there is intrinsically quenched disorder (the lists, see below), many such stable solutions are actually meta- 
stable in the sense that their cost function is not minimal. If the rules for breaking pairs are then softened, the society 
will go to a state of lower cost, i.e., more happiness. In some sense this is just what happens in reality. Indeed, due to 
changes in morality, marriage rules are weakened in the course of time. One century ago it was in most societies both 
morally and economically almost impossible to separate married couples, whereas nowadays it is a generally accepted 
£ — ' phenomenon in most Western and several Eastern societies. In line of the principle of freedom of the individual, many 
ON ■ American movies and television soap series nowadays even advocate that a person has the right to improve his/her 
own situation, no matter what cost this brings to the present partner, to the surroundings of the aimed new partner, 
or to society. To give a proper description of the total happiness of society, one should then, however, also take into 
account the cost for breaking up existing families. Such a more complete description would, naturally, put limits to 
. individual freedom. 

Let us consider a group A of N men and a group B of M woman. In reality the ratio v = M/N will slightly 
exceed unity. The man i has a wish-list, that is to say he evaluates all the women in group B and ranks them into 
a preference order. He will prefer, of course, the number one on his list, but that woman may not have reciprocal 
. intentions. Our task is to find a best global arrangement in which the total "happiness" is optimal. This calls for a 
cost function which we shall refer to as the Hamiltonian. In each arrangement preferably one man is paired to one 
woman. However, the possibility of remaining single is needed if the groups have a different number of members. Also 
polygamy can be permitted with assigned penalty weights. 

Globally stable solutions can be only found for small systems using exact enumeration methods. It is not hard to 
imagine that, as the system becomes larger, many conflicting wishes make a satisfactory solution harder to find. A 
systematical analysis is called for, which is conveniently done in terms of a partition sum at a fictitious temperature 
T; finally the T — ► limit should be taken. There is a connection with spin glasses where many different but 
roughly equivalent states occur. It was this connection which led in the mid-eighties to the replica analysis of several 
combinatorial problems Q , such as the bi-parted matching problem by Orland || , the mono-parted matching problem 
and the traveling salesman problem by Mezard and Parisi ||. 

The wish list can be implicitly given for the man i. For a given distribution he draws a (random) set of real 
positive numbers £fj, that gives the energy cost for pairing with woman j. This set can in principle be arranged into 
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a increasing order, with the most preferred partner having the least cost. Likewise, each woman j has a different 
(random) list giving the energy cost to form a pair with man i. 

In standard studies the cost of a configuration of marriages is the sum of rankings on the lists. For matching two 
groups of N members each person therefore brings an integer cost C, 1 < C < N, that equals the rank of his actual 
partner on his wish list. These integer costs are not so realistic, however. Number ten on a large list of candidates 
will almost be just as acceptable as number one, and not ten times less. For this reason we will follow Orland and 
Mezard-Parisi, who have considered non-integer costs, iy. 

The total cost for pair (i,j) is then £y = atfj + (1 — a)t\j. Note that in general £° ^ since intentions are not in 
general reciprocal, a characterizes another asymmetry of the problem: when a = 1 the cost function only represents 
men's interest, with that of women totally ignored, vice versa a = women's wish prevails. Besides these extremal 
men or women dominated situations, the case a = 1/2 represents the symmetrical (ideal) balance between the two 
groups. In society a is slightly tilt to men's advantage. 

Apart from this, if the man i remains single, this costs an energy A4. Basically, if Ai > mim, £fj then it is profitable 
for him to form a pair. Likewise, the cost for woman j to remain single is Bj. We shall also allow for polygamy in a 
grand-canonical fashion. If man i marries k > 2 women, there is an additional energy cost (k — where nu plays 

the role of a chemical potential. Similarly there is a cost (k — l)/i2j when woman j marries k > 2 man. Also the /it's 
can be taken as random variables, or may just be fixed. 

The goal of the problem is to find the best possible solution, i.e., the one that has lowest energy ("cost") for the 
whole society. It remains to be seen in how far this is the best solution for a typical individual, and which dynamics 
has the best state as stable solution. 

For equally large groups (M — N) and and no singles {Ai = B = oo) nor polygamy (fj,u = ji2j = oo) this marriage 
problem is closely related to the bipartite matching problem studied by Orland |f| and by Mezard and Parisi || . The 
major difference lies in the asymmetry of energy costs 1% ^ t\j, the new features such like being single or polygamy 
presents also new questions. Our analysis will follow these approaches as much as it is possible. 

The setup of the paper is as follows. In section II we describe some simple scaling aspects of the problem. In section 
III we describe a detailed replica analysis. In section IV we make a low T analysis. We close with a discussion and 
summary. In the appendix we consider the energy fluctuations. 



II. SCALING ANALYSIS OF THE PROBLEM 

A state is fully characterized by the set of numbers {riij}, where n,j = 1 if the pair is formed and zero else. 
The total number of partners of man i is rij = Y]j riij. Its typical value is rii = 1 (man i finds a woman/a woman 
finds him); however, rii = (man i remains single) or m > 2 (he is a polygamist) may also occur. In a similar way 
uij = J2i n ij describes the fate of woman j. The interesting case is where n, and rrij are annealed variables, so that 
both possibilities can be weighted thermally. 

The energy of a state is 

oo oo 
ij i i k—2 j j k—2 

Let us first consider the scaling for the pairs formed. The cost for pair (i,j) is 

^.=a4 Q) + (l-a)4 b) (2) 

In order to form a pair, this combined value should be as low as possible. First consider a = 1. Then typically the 
lowest i\y is £^l n is given by 
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PaWdlK- (3) 



Since we shall take 



Pa(i) = " , u (4) 

T(r a + 1) 



this implies = [r(r a + 1)/N] 1 ^ r ' l+1 \ The pairing energy would be equal 
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CW. = M^~JV^ /(r ' +1) (5) 

Typically, a man would find a woman that is high on the top of his list. The women would accept this, since their 
wishes have no effect what so ever. 

Now, if woman's preference is also taken into account (a < 1) things change. For a = 1/2 the distribution for the 
total I is: 

P(f) = r(r + 1) r = r a + r b + 1 (6) 
while for general a the small-^ behavior is 

^ M or. + i(i-a^r(r + i) ( ^ 0) (7) 

The typical value £t yp ~ tin ~ VlH-i) now implies that the typical partner is located on position k of the wish 
list related as ^ tyF d£p a (£) = k/N, yielding k ~ N'( r t>+ 1 )/( r + 1 ) . For the symmetric case r a = r& = (r — l)/2 this 
becomes fc ~ vW. So for the two-side weighted case the wishes of each individual are quite moderately fulfilled: the 
ideal partner is not within reach. She/he has typically has different wishes, which excludes formation of the pair. For 
the bulk of pairs the partners are only moderately biased towards each other. This is the price society has to pay for 
individual freedom. 

The ground state energy will scale as Uq ~ Ni typ ~ N r /( r+1 ). This implies that bachelors and polygamist s only 
play an interesting role when their costs (Ai, Bj, fi\ t i and fJ>2,j) scale as Uq/N ~ iV — jf ^ ne cos ts have a larger 
magnitude, these options do not occur provided the groups have equal size; if they are smaller, no pairs will be formed. 
Notice that these scalings are not universal to the individual, but depend on the exponent r entering the distribution 
(fjj) of costs for pairings in the whole societies. 

III. THE PARTITION SUM AND REPLICA'S 

We consider the problem in a statistical mechanics approach. We thus consider the sum over all configurations 
of Boltzmann factors exp(— (3E) at given temperature T = 1/(3. The optimal solution has lowest energy; it will be 
dominant in the limit T — > 0. The statistical approach also yields the entropy, related to the degeneracy of the optimal 
solution, as well as the average numbers of married man etc. 

To simplify notation we first take sure values for Ai, Bj, fin, fi2j- The partition sum can be written as 

N M 1 oo 

^nnE e 

i—l j — 1 riij—0 m ,rn,j—0 

1 00 M 

= V" e^£r=i {(A+w)i "*' 0+(ni " 1) '' l} e"^i=> {(B+, ' 2) ' i °' j,0+(mrl) '' j} 

riij—0 ni.rrtj—0 

N M M N 

e _/3E « niilii II 6i "- - E II 5 ( m i - E "«) ( g ) 
i=i j=i j=i i=i 

where the 6's are Kronecker ^-functions. Repeatedly using the integral representation 

6(a-b)= r^e^'-Q (9) 
Jo 27r 

one arrives at a form where the sums can be carried out. After doing that one gets 



z = TT / lit } —T\ / — - '1— TT n +e~ ,A * _ w~^l do) 

i=l 3=1 *,3=1 

In order to calculate the quenched average free energy one employs the replica trick: 
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F = -Tin Z = -T lim 



Z n - 1 



n—>0 n 

One thus replicates the partition sum n times and uses for each pair 



(11) 



no 



') = i+E E 



-i(A° 1 +--- + A° p )-i(M° 1 +"-+' 1 'i >P )-j ) << 3 ^. 



£*=1 



1 l<ai<---<a p <n 



Now the quenched average over £y can be carried out. One gets 



= / d^.p(^.)<p(4) e " P/S( ° 4&+(1 " 0,) ^ ) = | 



This defines 
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(1 + ap)8) r »+ 1 (l + (1 - a>/3)^ +1 (p^^a^+^l - a)' r " +1 



where r = r a + rj, + 1 . In terms of r our g p has a similar form as in the matching problem || Q . 
We shall confine ourselves to low temperatures, where we scale temperature with N such that 



T = T 



N 



l/(r+l) 



is of order unity. This implies that 
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r+l 



is also of order unity. We then get to leading order in N 

Z*-TT ( dX " e ~ PA + £f ~ e~ /?(A+A ' l) )e' A " -pr /• e-^ B + (1 - er^+^)y^ 



2?r 



n 



2tt 



1 - e _/3/ia+i/1 " 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



with 



1 - 

5 = ^EE^ E 



3 _ 1 (a°i+... + a° ! ')-«(m7 1 +-+^ ! ') 



*J P=l l<ai<---<a p <n 

The i, ^-dependence of expS can be decoupled by a Hubbard-Stratonovich-type transformation 



n „ 

pS = II 11 / dP ai ... ap dQ ai ... ap exp 



P=l {«r} 



p=l ^ p l< ai < - <a p <n 



N 



M 



ex?E E p E e^ (A " 1+ - +A ° P) + Q Q ,.. Qp E e ~ 

p— 1 l<Qi<---<a p <n z— 1 j—1 



t(/i° 1 +-+/i° p )l 



(18) 



(19) 



Now all i,j indices are decoupled. The Xi integrals all yield the same factor. The replicated partition sum thus 
becomes equal to a P, Q integral of a function exp[— N$(P, Q)]. Because N is large, we may approximate the integral 
by its saddle point value. Therefore the replicated free energy F n = —TlnZ n is just equal to this saddle point value. 
We thus obtain 



^-V— V O P 



N i ' Or, ^ 

p=l JF l<cti< - <a p <ra 

where the latter two objects are defined as 



(20) 
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and similarly 
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n 
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dX a e 



-PA . 



(1 



-0(A+iii)\ e i\ 



2tt 



I — e —0tMl+i\ o 



/ j ■ r ai—a p c 
p—l l<Qi<-"<a p <n 



+ •••+>„,) 



(21) 



2tt 



d^ a e-? B + (1 - e -^(B+M»)) e i/*a 



2tt 



1 — g— /9A*2+»/i a 



<E E 



1 l<CKi<---<a p <ra 



-ifMojH hMo p ) 



(22) 



Following Orland and Mezard-Parisi we assume that the relevant saddle point has replica symmetry, viz. Q ai —a P = 
Qp, Pon - otp = Pp- Then the X a integrals can be replaced by contour integrals over z a = exp(— iX a ), and the integral 
can be evaluated from the poles at z a = and exp(— /3fi±). This yields 



- f>A - f> " 1 )8*c,,o + S, a ,e-»i\expY, P P E 



p—l l<cti<---<a p <n 
-0H1P 



= e"^ 1 E (fc) (-1 + e- /3A - /5AI1 ) fe exp^ P p 

/c=0 p=l 

e* £ Q (-DM - 



0) 



-Pmp 



T(k -n- i0)p\ 



(23) 



where iO = it with e — * first. In the last step extension of the sums to oo was allowed because the extra terms are 
identically zero. So far n has been integer. Now we can continue this expression for n — > 0. There are two types of 
contributions: the k = term equals 1 + 0(n), while the k > 1 are all 0(n). Indeed, for those terms one uses 



T{n + i0+l) T(n + l)sin(Tr(k-n-iO))T(k-n-iO) n(-l) 



fe-i 



T(n-k + l)k\ 
In the limit n — > one thus gets 



-(l-nV>(fc)+n^(l) + ---) (24) 



z^^nz a +n 2 z^ +0(n 3 ) 



(25) 



The result for z a is presented in eq. (|28|), while is presented in eq. (|7< 
Likewise z& = lim„^o z n/ n ^ given in eq. (^9|). 
Further we find, using eq. (|2~i|), 



QpPp ( 



p =i 3p i< Ql <... <Qp < n p =i 3p VP 



<E 



(-l)P- 1 Q P P P 
P 9p 



(26) 



A. The quenched free energy 

Now the limit n — ► can be taken simply. The quenched average free energy F = lim„^o F n /n thus follows from 
eq. © as 

/3F ^ Q p P p (-IF' 1 

-TT = 2^^^ z a(AMi) - PZ b {B,fj, 2 ) (27) 

JV p =i 5 p p 

with 

^, w) . + g (zD^-^. _ £ g - ^"')' exp {g Pp (zlg|±£) e - rfm } (28) 

and 
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p=i ^ fc=i ' U=i v J 

In section 2 we have discussed the possibility that A», Bj, /ii^ and are independent random variables. For that 
case the analysis goes along the same steps; one only has to average the present expressions for z a and Zb with respect 
to these variables. 

In the case where these variables are non-random, and only the pair costs are random, the saddle point 

equations read: 



fc=l 



and 



where —P(k) and —Q(k) are the arguments of the exponent in the fc'th term of z a and Zj,, respectively. 

We can compare this with the Orland and Mezard-Parisi results. In the limit fii t 2 — > oo only large-fc terms are 
relevant. Introducing £ = ke~ P^ 1 - 2 we can go to an integral representation. Adding and subtracting to z a the 
expression X^fcLi ex P( — Q/k we regularize the small £ behavior. We find that the contributions cancel, and obtain 
in the limit 

* ft- P io ^ P p tt ? ! J 

yl e -exp 2^ -j e 5 ] (32) 

In absence of bachelors (A, B — ► oo) and the limit of equal group sizes (y — M/N — > 1) this equation reduces to the 
result of Orland and Mezard-Parisi. Equation ( |32"|) was also derived by us using their method of taking derivatives 
rather than summing over residues. 

The saddle point equations follow from eqs. (^), (jf^), and (p9|). 

Q, = tf^-^>-exp£^ (33) 

and 



P< Jo £=J Pi 



Pp = A r dte-^e- 1 expV ^Hl! (34 ) 
P- Jo p[ p! 



IV. LOW TEMPERATURES 



From eq. (15) one sees that g p = (p/3) ( 1+r ) is of order unity if we scale temperatures with N such that f ~ 
2 n 7V 1 /( 1 + r ) remains fixed. As discussed in section II, the free energy is not extensive but scales as F — —T\ogZ = 
^N-i/(r+i)f \ ogZ = AT/(r+i)^ with intensive p. 

For having non-trivial behavior one also needs that A scales properly with N, A = aN -1 '^ 1 ' so that (3 A = /3a, 
and similarly = fjli,iN~ 1 ^ r+1 ' . 

Inserting £ = exp((3£) we can now introduce the functions P(£) — P(£), etc. To simplify notation, we shall however, 
again denote P{1) by P(£), etc. Thus 

P m = y (- 1 ) p " li P e P^ QW = y (-i) y ~ x Qp e ^ (35) 
£1 p [ t^l p ] 
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At finite T they satisfy 

/oc 
dyB(£ + y)e- Q( ^e- e (36) 
-oo 

which for T — > reduces to 

P{£) = v J dyB(£ + y)e- Qiv) (37) 

Similarly 

rOO r-b 

Q{£) = dyB(l + y)e- p Me- e0(v ~ b) -> / dyB^ + y) e - p ^ (38) 



with 

(_l)p-i e J>W-2/) 



P =i 


.iy-i e ppe 
p r p\p\ 




/•OO 

/ dyj/V 
Jo 


r(r + i) 




£>0 


_^ I r(r+i) 
1° 


Bessel function. 





r ( r + 1 )^o ^ ^ p!(p-l)! 



(39) 



A. Observables 

The internal energy scales with TV as 

^=|^=l(a'-+ 1 (l-«r +1 ) V(r+1) ^ (40) 
The scaled internal energy £/ = d((3F)/d(3 consists of three terms 

U=Uab + Ua + Ub (41) 

The pair energy is 

u ab = T ^ 



/OO 
d£P{£)e^ m ~ a) e~ p ^ ^ (1 + r) / d£P(£) e - p ^ 
-OO J —oo 

/OO 
d£Q(£)e~ e e ^ -» (l + r)i/ / d£Q(£)e~^ (42) 
-OO J —OO 

The male-bachelors energy is 

j,Ua = a d£PeW- c h- p ^- eW - c) -+ ae- p ^ (43) 

■'" J —oo 

and the female-bachelors energy is 

l?7 b = 1/6 / d^^-^e-^^)-^^ 6 ' - i/fcr*< 6 > (44) 
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The average fraction of married man is obtained as 



= 1 - / d^V^^V^ -> 1 - e- # (°) (45) 



and likewise the fraction of married woman 



< m > M 

3=1 



l^( mj .)^l- e -«3W (46) 
j'=i 

The case of Orland, Mezard and Parisi is recovered for v = 1 and a = b = oo (infinite energy cost to remain single). 

B. A special case at T = 

Let us follow Orland and Mezard-Parisi and take r — 0. We assume equal costs a = b for males and females to 
remain single. (This restriction is not necessary; our results can be generalized to the case a ^ b.) For this case eqs. 
become 

dy e -P(.v) p(£) = v I dy e -QW (47) 
-i J —i 

So that P = Q = for I < -a. One gets 

P'(£) = „ e -«(-fl Q'(l) = e - p ^ (48) 
These equations can be combined in a single equation for P 

P"(l)=P'{l) e - p{l) (49) 
which can be integrated once. The solution reads 

P{1) = log H > ___ QW = i og[ _( + i)] (50 ) 

7 7 7 — J- 

where 7 is a parameter. It is obvious that P{—a) = 0, while Q{—a) — requires that 7 is the largest root of 

e- 27a = (7-l)(--l) (51) 
v 

From this we derive 

(n) = i/(m) = 1 + 1> - 7 (52) 
In the present case (r = 0) the energy U ~ O(N ) is non-extensive. The contributions to the energy are at T = 0: 

00 /„A-ft 



7 , /-, v,kMl-v) 1 



, k 2 7 fe fc 2 

fc=l x ' 



7 ln^ln( 7 -^)+ 7 f; 1 (53) 



Lu a = a( 7 -v) 1^ = 0(7-1) (54) 



We can now check some simple cases. 
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• For a — > oo and v = 1 one has the Orland-Mezard-Parisi case with 7 = 1, so all men and women are coupled. 
If we choose a = 1/2, we recover the result U a b = 7r 2 /12, while there is no bachelors energy. 

• If a 3> 1 at fixed v > 1, it follows that 7 = f, yielding and (n) = z/(m) = 1 and C/ a = 0, saying indeed that all 
TV men are coupled to N out of the M — vN women, a fraction \/v of the total. There is no way to avoid the 
energy cost £/& = a(v — 1) = (M — N)A of the unpairable M — N women, where A is the energy to pay for each 
bachelor. Likewise, for v < 1 there are more men than women, and the result 7 = 1, (l/v)(n) = (m) = 1 says 
that now all women are coupled to a fraction v of the men. This time the surplus of males is responsible for the 
energy cost (N — M)A. 

• For small a it holds that U ~ U a b ~ (n) ~ a. In the appendix we shall find that the fluctuations scale as 
SU ~ a 3 / 2 /\/N, implying 5U/U ~ y/a/N. As expected, it vanishes for a — > 0. 

• In the limit a — * it is energetically more advantageous to remain single. Though pairing may occur at finite 
T due to entropic reasons, the T = result 7 = 1 + ^ indeed implies (n) = (m) = 0, U a t, = 0, expressing that 
no pairs are formed at T = 0. 

• Let us finally consider equally large groups (y — 1) and finite scaled cost for remaining single, a = b = 0(1). 
Then it follows that 7 = l + exp(— 7a). The result < n >—< m >= 1 — exp(— 7a) expresses that a finite fraction 
of bachelors remains, because it is more advantageous that this fraction of the individuals is not paired. This 
aspect remains in the more general case v 7^ 1 and a =/= b. 

1. The entropy at low temperatures 

A test for checking the correctness of the replica symmetric approach is to calculate the leading behavior of the 
entropy at low temperatures. For doing that, an expansion in powers of T is needed. 

For simplicity we consider the case r = 0, a = b = 00, v = 1. The kernel B from eq. ( |39| ) can be expanded in 
distribution sense 

B{t) = 1 - J (2e/ M/2 ) = 6{£) + ai fS(£) + a 2 f 2 5'{£) + a 3 f 3 S"(£) ■ ■ ■ (55) 
Partial integrations reveal that 

/•oo />oo 

ax = I d£ {B(£) - 6{£)) = -2 / da; Ji (a;) log ~ (56) 







Similarly 



and 



00 _ 

x. 



a 2 = -P z I d££{B(e)-6(£)-axT5(e)) =2 / da; Ji(x)(log -) z (57) 







a 3 = -/3 3 / d££ 2 (B(£) - 9(£) - ai T8(£) - a 2 5'{£)) = -- dx Ji(a;)(log -) 3 (58) 
1 J -00 ^ Jo * 



It holds that 



Jo r(i-| M ) 



(59) 



Expansion in powers of /1 yields the needed integrals. One obtains 

a!=27 B ; a 2 = 2 7 |; a 3 = ^(3) + 1 7 | ( 60 ) 

where 7^ = 0.577215 is Euler's constant, and £(3) = X*Li V^ 3 i s a Rieman zeta function. 
Eq. (pq) reads for the present case 
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/oo 
dyB{£ + y)e- p ^ (61) 
- OO 

Expanding 

P(y) = P (y) + TP l {y) + f 2 P 2 (y) + f 3 P 3 (y) + ■■■ (62) 



we find using ( p5| 



Pi(£) = - I dye- p «^P 1 (y) + 2 lE 



-i 



/CO -i 
&ye- p ^y\~P 2 {y) + -Fi(y) 2 ) - 2 lE e- p ^ + 2 7 % e - p °^ P(,(-£) (63) 



The solution reads 



P (£)=log(e« + l) (64) 

Pi(i)=lE^—r (65) 
e l + 1 

p « = ^(?TiF (66) 



From eq. we obtain the entropy 



= pj°°d£ \[9{-i) - e- p °^] + [e~ e$t - 9{-£) +f e - p «^P 1 {£)}+f 2 [e- p °^(P 2 {£) - ^Pi(£) 2 )]^j (67) 

The terms calculated so far lead to orders T , T° and T. It turns out that all three prefactors vanish. Thus S is at 
least of order T 2 . We have to go one step further. The equation for P 3 reads 



P 3 {£) = / dye 



-Po(y) 



• + y){-Ps(y) + MvWv) - l p i(v)} + 2 ie8(£ + y){-P 2 {y) + \p?(v)} 
+2 7 %S'(£ + y){-Pi(y)} + (|c(3) + \l%)5"{£ + y) 



(68) 



Some tedious analysis reveals that the solution reads 



e 



P3 W = -T ?TT " f FTTF + (C(3) + i 7l) (?TTF (69) 



The next order contribution to the entropy reads 



It follows that 



O poo 1 

^=T 2 _ dye- p °W{P 3 (y)-P 2 (y)P 1 (y) + -P?(y)} = -T 2 P 3 (ac) (70) 



S = \n((3)T 2 (71) 

The specific heat therefore equals 

C = ^C(3)T 2 (72) 

Notice that in S all contributions of order 7^ have canceled. Therefore we could have simplified the calculation by 
neglecting je- This makes it possible to extend the calculations to the case where a is finite. 
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C. Random costs for bachelors and T = 



If the costs for bachelors are random with scaled densities pi(a) and p 2 {b), then eqs. ( E^ ) become 

dyiix(y)e- p W P(£)=v R 2 (y)dye-^ (73) 

where 

/■oo /-oo 

= / da Pl (a); R 2 (y) = / dap 2 (6) (74) 



This leads to the differential form 

Q'(£) = i? 1 (-£) e - p (-^; P'{1) = V R 2 (~£)e- Q ^ (75) 

and 

P"(i) = -J^f/'W + Ri(£)e- m P'(£) (76) 

We have not been able to solve these equations for broad distributions p\ and p 2 . In case that the costs Ai and Bj 
for remaining bachelor take discrete values only, pi(a) and p 2 (b) are sums of delta- functions. Then the solution can 
be constructed by generalizing our previous solution in the segments between the (5-functions. 



V. DISCUSSION 



We have considered the bi-partite matching problem, where members of each group have a wish-list of non-negative 
real numbers ranking the possible partners from the other group. This is a polynomial (P) problem, for which fast 
algorithms are available. We have extended previous approaches || [|| to the case where groups have different size, 
members can remain single, and polygamy occurs. Our results can be extended to the case of random costs for 
remaining single or being polygamic, by averaging eq. ( p7j ) over these variables. 

In section 2 we have given a scaling analysis for the problem. The interesting regime is where the £'s, the costs for 
pairings, are low. We assume that the probability density of the costs of males ft a \ scales as a powerlaw, (£^) ra , 
and the one of females as (£^) r ». It is a standard result of statistics that the probability density for finding a couple 
with low cost £ = &^ + £( b > then scales as £ r with r — r a + + 1. The exponent relation r + 1 = (r a + 1) + {r^ + 1) 
expresses that both partners must have a low cost. 

The typical cost can be estimated from = 1/iV. Therefore the energy will scale as U ~ N£t yp ~ ]\f r /( r + 1 ). 
Bachelors and polygamists only play an interesting role when their costs (Ai, Bj, and p 2 j) scale as U/N ~ 
7Y- 1 /(»+ 1 ) jf these costs have a larger magnitude, such options do not occur (provided the groups have equal size); 
if they are smaller, no pairs will be formed. Notice that these scalings of costs for bachelors and polygamists are 
not universal to the individual, but depend on the exponent r entering the distribution £ r of costs for pairings in the 
whole society. 

We have not worked out the role of polygamy at low temperatures. As it must play a role in case of unequal group 
sizes (v 7^ 1) but absence of bachelors (A = B = 00), we expect it generally to play a non-trivial role down to T = 0. 

In a special case at T = we generalized previously known exact solutions to include for bachelors and unequal 
group sizes. It could be checked that if the typical cost for being bachelor is of the right order of magnitude, quite a 
number of members of the groups will remain bachelor, even if the groups have equal sizes. 

It is found that the entropy remains positive and vanishes quadratically in T. This result supports the validity of 
the replica symmetric approach. 

A calculation of the variance of the free energy is performed in the Appendix. In principle replica symmetric 
approaches may lead to a negative (and thus wrong) prediction for this variance. The somewhat painful analysis 
shows that the replica symmetric Ansatz produces the correct sign. In the totally symmetric case of equal group sizes, 
equal weighting of males' and females' wishes, and in the absence of bachelors or polygamy (with parameters v = 1, 
a = 1/2, A = B = p x = n 2 = 00, r a = r b = -1/2) we derive SU 2 = 0.442878/7V. It has the expected scaling in 1/N. 
The correctness of the sign supports again the expectation that replica symmetry yields the correct result for this 
problem. In its turn replica symmetry is expected to lead at worst to a polynomial degeneracy of the ground state. 

Our results can be summarized by saying that whether bachelors occur or not, depends not only on their intrinsic 
capacities, or lack thereof, but also on global aspects of pair formation in the society. 
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A. APPENDIX: VARIANCE OF THE FREE ENERGY 



The present problem has only been considered within the replica symmetric Ansatz. To test the validity of this 
approach several checks are possible: a stability analysis, the sign of the entropy or the sign of the second cumulant 
of the free energy. 

The physical free energy should have a positive variance. This implies that the leading term in N should be non- 
negative (if the leading term vanishes, then the fluctuations have typical amplitude less than \/~N). It was pointed 
out by Saakian and Nieuwenhuizen, that the 0(n 2 ) of F n can be used to check non- negativity of the variance M. 
Indeed, for F = -TlogZ and F n = -TlogZ", it holds that 

F n - n~F 



f3 (F - F) 2 = -2/3 lim -2— (77) 

n^O n 

In their variational analysis of an interface in a disordered medium, Saakian and Nieuwenhuizen found that the 
variance is negative if no replica symmetry breaking is taken into account. For a finite number of breakings the 
variance becomes smaller in magnitude, but remains negative. Only for an infinity of breakings it becomes zero to 
leading order. 

For simplicity we consider here the symmetric case A = B, v = 1, /ii = /i2 = ft- We first need the n 2 term of z®, 
Za , defined by eq. (p5|), 



(-V)PQ e -Pw e-QWd - e -P(A+mh* 

zf) = £ — {m _ ^(i)) + £ £ )_ {m _ 



~ e -OW(l _ e -/3(A+ Ml ))fc ~ ( _ 1)P -iQ pe -ff w i /~ e -Q(k)( 1 _ e -p(A+ m )- } k \ 2 

^ k ^ p 2 I ^ k j 

fc=i p=i r \fc=i / 

E - — { — k L E i-ir-^e-^^^mk+p) - m) m 

k=l p=l ^ ' 



From eqs. @, ©, ©, ©, @, and © the variance follows as 



«2 i 00 (_-[\PQ2 

JW 1 ^ = - \ E ^^WP) - *(D) + *£° (79) 



In the first term we insert the equation of motion (p0[). We consider the obtained expression in the limit /i — > oo. As 
before, the only possible surviving terms are those where e - ^ can pick up a factor k. The relevant terms are 

«2 i °L (_i \pq p-pVti-nfi. + n \ P -Q{k)(i _ e -p(A+n)\k 

P_ iF2 _ _ i £ j i)o^r ( . + P ) w(p) _ „ (1)) !^_(L_f L 

+ g, e^Hi - «-"'«-)' wt) _ _ i f£ «-«*a - V (80) 

fe=l \fc=l / 

In the limit /x — > oo one obtains an integral representation involving the variable £ = fce - ^. We again add and 
subtract to e~^^(l — e~^(' 4+A11 )) fe to regularize the small ^-behavior. This leads to 



iN 2 Jo 6 P ! 7o £ 

+ - 7^) + - + [ f(e-^e-^ A - e~*)\ % (81) 
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where the last Euler's constant arises from —-0(1), while the first comes from an integral representation for the 
sum involving ip(k). Notice that all /x-dependent terms cancel, confirming that the situation without polygamy is well 
defined. Using one has finally 



4/V" 



SF 2 



e -QW e -te-* A 



drj 
o V 



2 Jo € 



where 



^^(^)f(^(p)-^(i)) 



(82) 



(83) 



As can be seen from the first non-vanishing contribution, the p = 2 term, £?2 is a positive function. 

We now go to T — 0, so the internal energy equals the free energy. We calculate B 2 for the case r = 0. We have 
the scaling T ~ T/JV, A = a/N. We set £ = exp/3^, 77 = exp/3f and insert in Q the identity 



V>(p) - ^(l) 



dz ln(l — z) 



2?rizP 1 - z 

This brings us to a form where the p-sum can be carried out, yielding 



B 2 {l + £')=f 



dz ln(l - z) 



1- Jo 



(84) 



(85) 



2iri 1 — z 

We can deform the z-contour and remain with the discontinuity along the real axis z > 1, which yields for £ + £' > 



B 2 (£ + £')=f 



/■°° dz 


/ 2e /§(W)/2y 


/a 2-1 




r Jo l ^ JJ 



■O(f) 



£ + £' + 0(f) 



(86) 



while B 2 = for £ + 1' < 0. 

We calculate the variance of the internal energy for the simplest, solvable case r — 0, a = b. To leading order, all 
terms in (|82|) are proportional to jj 2 ~ f3 2 /N 2 , implying 



SU 2 = {-h + 2I 2 - II) 



(87) 



Since U ~ N r ^ r+1 ^ = N° it follows that SU/U ~ N~V 2 , as usual. It holds that 



h 



d£ / d£'e-^-^ (£ + £') 



and 



and 



d££(e 



-Q(i) 



1)+ / d££e'^ 



d£{e 



-QW 



1) + / dte~ w) = (7 - l)a 



(89) 



(90) 



For a — > 00 (no bachelors) one has 

f°° 1"°° IT 2 

h= de\n(l + e e ) ln(l + e- f ) = 2* 1.202056; h = 2 \ lnfl + e~ l ) = — = 1.644934; 7 3 = (91) 
J-oo Jo 6 
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Since r = 0, boundedness of the entropy at low T implies that in this limit ^3 = 0, as can also be checked explicitly. 
Therefore the combination —I\ + 2I2 — Z| is positive. In the totally symmetric case (equal group sizes; equal role of 
males and females) with v = 1 a = 1/2, r a = rb = —1/2 we thus find 

(92, 

(Note that this implies for mono-parted matching problem with r = 0, p(£) = exp(— £), considered by Mezard and 
Parisi f|: SIP = 1.771512/iV). For a -> it holds that 7 = 1 + e~ 7a - 2 - 2a, and the prefactor 

- h + 27 2 - 1% + 2(y - la 3 ) -(a- 2a 2 ) 2 = ^ + G(a 4 ) (93) 

is still positive. 

It is thus seen that the replica symmetric prediction for the (free) energy fluctuations has the correct sign. This 
supports the expectation that our replica symmetric solution is exact. 
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